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Spin-dependent ballistic transport through a tunnel barrier is treated within the 
one-dimensional nearly-free-electron model. The comparison with free electrons 
reveals significant effects of band gaps, in particular in the bias dependence. The 
results are qualitatively explained by the number of incident and transmitted 
states in the leads. With an extension to ferromagnetic leads the bias dependence 
of tunnel magneto-resistance is discussed. 

Keywords: Nearly free electrons; Ttinneling; Ballistic transport; Bias dependence; 
Tunnel magneto-resistance; 
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1. Introduction 



Spin-dependent transport of electrons in nanostructures is one of most rapidly evolving 
areas of contemporary physics [1,2]. Important contributions to the field come from ap- 
plied physics, focusing on design and characterization of devices and applications, and 
from experimental physics, aiming at the understanding of fundamental effects. Contribu- 
tions from theoretical physics comprise on one hand transport calculations that are based 
on state-of-the-art electronic-structure calculations. By this means, properties of specific 
systems are investigated. On the other hand, simple — and sometimes too simple — model 
calculations are performed in order to explain experimental results, thereby neglecting 
occasionally important aspects of the system under consideration. Apparently, there is a 
considerably large gap to be filled between the advanced electronic-structure and the model 
calculations. 

The giant magneto-resistance (GMR) had and still has a huge impact on industrial 
applications [3]. The tunnel magneto-resistance (TMR), however, is expected to have 
even more potential for future applications and devices [4,5]. From a theoretical point 
of view, TMR has the advantage that it can be understood in terms of simple quantum- 
mechanical effects (e.g., the tunnel effect) and can be modeled rather easily. Therefore, it 
lends itself support for an investigation which connects model calculations to the frontiers 
of nanoscience. 

The purpose of the present paper is to provide a theory of spin-dependent transport 
through tunnel junctions which bridges the aforementioned gap. Its framework is the 
Landauer-Biittiker theory [6], in which elastic and ballistic transport is viewed as trans- 
mission of scattering channels through the device. All ingredients needed for the calcu- 
lations of the tunnel conductance can be computed step-by-step, therefore allowing for 
simplifications or extensions. The theory outlined here was implemented one-to-one in a 
set of Mathematica"'"' notebooks [7], in order to obtain numerical results. Hence, the 
properties of the tunnel junction can be easily manipulated, results for different set-ups 
computed and visualized rapidly [8]. 

Model calculations typically assume free electrons in the metallic leads that are con- 
nected to the tunnel barrier, whereas electronic-structure calculations deal with a much 
more complicated bandstructure. An important feature of the latter — which is missing 
in the former — are band gaps, i.e., energetic regions in which electronic states are not 
present. Obviously, these have a significant impact on the electronic transport. In order to 
take band gaps into account, we go beyond the free-electron approximation by applying the 
nearly-free-electron (NFE) model. Assuming one-dimensional leads, the theory is kept sim- 
ple enough to be tractable semi-analytically but it contains already the relevant aspects of 
the electronic structure. Both ferromagnetic and nonmagnetic leads are discussed, because 
the results for the latter can easily be transfered to the magnetic case. We note in passing 
that tunnel junctions with ferromagnetic leads typically comprise transition metals (Fe, 
Co, Ni), the electronic structure of which cannot be described well within the free-electron 
model. 

The paper is organized as follows. First, basic theories for ballistic transport which are 
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relevant to the current work are sketched (section ^ . Having introduced the model for 
the tunnel junction (section O]) , the electronic-structure calculations (section H]) and the 
treatment of tunneling (section E]) are presented. Representative results of the tunneling 
calculations are discussed in suhsection l6.2| focusing first on tunnel junctions with nonmag- 
netic leads (subsection 16. Ij] . Spin-dependent transport is investigated in subsection 16.21 

2. Basic theories for ballistic transport 

The theory presented in the forthcoming sections applies to elastic and ballistic transport 
of electrons through a tunnel junction (figure HJ. Hence, the method of wavefunction 
matching for calculating the scattering wavefunction in the entire tunnel junction can be 
applied. Experimentally, ballistic transport is observed for small and defect-free samples. 

The simplest approach to tunnel magneto-resistance traces back to the work of Julliere 
[9], where TMR is related to the polarization of the junction defined via spin-dependent 
numbers of states at the Fermi energies Ep in the leads. Maekawa and Gafvert [10] modified 
this model by defining the polarization in terms of the densities of states at Ep. 

In the Landauer-Biittiker theory [11], the conductance is proportional to the transmis- 
sion probability of the scattering channels ("conductance by transmission"), G oc T^Ep). 
For a tunnel junction, the scattering channels are the eigenstates of the leads in absence of 
the tunnel barrier, i.e., the Bloch states. This approach can be applied in both sophisti- 
cated electronic-structure calculations or in model calculations. For example, Slonczewski 
assumed free electrons in the ferromagnetic leads and a step-shaped tunnel barrier (as 
sketched in figure HJ [12]. 

The present work goes beyond Slonczewski's model, in that it takes into account the 
nonzero potential in the leads. Applying the NFE model, the occurrence of band gaps in 
the electronic structure has a profound effect on the electronic transport. For the tunnel 
barrier, we stick with the step shape, being aware that a differently shaped barrier will 
lead to quantitatively different conductances [13]. The qualitative picture, however, will 
still be valid [14]. Although being computed in the framework of the Landauer-Biittiker 
theory, the present results can also be interpreted in terms of the densities of states. 

3. Model for the tunnel junction 

The one-dimensional tunnel junction consists of three regions: the left lead £ = {x\x < 
-d/2}, the barrier B = {x\ - d/2 < x < d/2}, and the right lead 7^ = {x\d/2 < x}. To 
keep the model as simple as possible, the leads were chosen identical, i.e., identical lattice 
spacing and potential (This restriction can, however, be relaxed in the Mathematica™ 
notebooks [8]). Within the leads, atoms are positioned at x = —{d/2 + no) and x = 
+ {d/2 + no) for C and 71, respectively (with n semi-positive integer). The lattice spacing 
a serves as the unit of length. 

In accordance with the above spatial decomposition, the potential of the tunnel junction 
consists also of three parts (figureEj), namely V'-' , V^, and V^. The barrier potential V^{x) 
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is assumed constant, with value Vbarr > 0. In principle, any other shape could be used [15]. 
In the leads, the potential is periodic, i. e., V'-'{x) = V'-'{x — a) and V'^{x) = V^{x + a). 



4. Calculation of the electronic structures 



4.1 Electronic structure of the leads 

According to the Landauer-Biittiker theory, one needs to compute the scattering channels, 
i.e., the electronic states of the leads in absence of the tunnel barrier. Thus, the lead 
potential extends over the whole x-axis. 

The Schrodinger equation Hip = Eip in atomic units [16] reads 

i){x)=E^{x). (1) 

Expanding the periodic potential into a Fourier series, 

\/(x) = ^V?^ (2) 

geQ 

the requirement of periodicity \V{x) = V{x + no), n integer] restricts q immediately to 
the reciprocal lattice Q, i.e., g G Q = {ng\n integer, = 27r/a}. The mean value Vq of the 
potential serves as origin of the energy scale and, hence, is set to zero (Vq = 0). Assuming 
further a real V{x) and the mirror symmetry V{x) = V{—x), one has V^g = Vg Wq E Q. 
For free electrons, V{x) vanishes identically {Vg = Vg G Q). 

For a periodic potential, the eigenstates ip of the Hamiltonian H fulfill Floquet's theorem 
[17]. In analogy to (equation EJ, the Bloch states are as well expanded into a Fourier series, 

^,(x)=e^'=^-^c,(A;)e^'^^ (3) 

In the following, the Schrodinger equation is solved by two methods, which have in 
common that the differential equation is transformed into an algebraic set. 



1 d^ 



4.1.1 Computation of E{k). Inserting equation ((2j and equation ^ into the Schrodinger 
equation ((TJ gives the infinite algebraic set 



,ik) + J2Vg-g'Cgik) = Eik)cgik), Vge 
g'eQ 



(4) 



with k e]- g/2,g/2]. Defining 



Hgg'ik) = ^iq + kf5gg, + Vg^g,, V^, q' e 



(5) 
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one obtains the compact form 

J2 H,,.{k)c,,{k) = E{k)c,{k), Vg G Q, (6) 

which can be cast into matrix form, H{k)c{k) = E{k)c{k). The Hamilton matrix H{k) is 
square and hermitian, resulting in real eigenvalues E{k). 

4.1.2 Computation of k{E). In tunneling calculations, an electron is specified by its 
energy E. Since the Bloch states in each lead are characterized by the wavenumber k'~ and 
A;^, respectively, it has to be assured that E{k'~) = E{k^) = E. Alternatively, one could 
solve for k{E) in each lead, an obviously favorable method [18]. 
Rewriting equation (jH) as 

(^q -^J2{E + Va) + k^[q+ ^2 {E + ^) + A:) c, + 2 ^i-i'^'i' = 0' ^9 e Q, (7) 

an algebraic set of double rank as in the E{k) approach is obtained by defining auxiliary 
expansion coefficients Wq-. 

[q -^2{E + Vo) + k)wg + 2 J2 V<i-g'Cq' = 0, 

9'6Qxg (8) 

-Wg + (g + {E + Vo) + k)cg = 0. 

In matrix form, equation ^ reads 

hn{E) hu{E) \ ( w{E) \ _ ( w{E) \ 
h^,{E) h22iE) J [ ciE) J -^^''^ \ c{E) ) ' 

where the Hamilton matrix H{E) is a 2 x 2 supermatrix comprising submatrices hij of 
rank ro, 



= {q- ^2 (E + Vo)) Sgg>, Vg, g' G Q, 



qq 



hf^^^ = {q+V'^{E + Vo))5,q,, 

"'qq' - -^qq'- 



(10) 



Since H(E) is not hermitian, its 2rQ eigenvalues k^^\E) can be complex. The dis- 
persion relation k{E) defines then the complex band structure [19]. For the tunneling, 
real fc*^*-' are appropriate because only these belong to square-integrable wavefunctions 
[t/'W G L2(— oo, +oo) for k^"^^ real]. If, however, Im/c*^*) > (ImA;*^*) < 0), the eigenfunction 
i/j^^^ decays exponentially towards +x {—x) [20]. 
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5. Calculation of tunneling 

Having specified the electronic states in the three regions £, B, and TZ, the wavefunction 
in the entire tunnel junction at a given energy E has to be determined. This is achieved 
by requiring continuity of the total wavefunction and of its first spatial derivative at the 
two interfaces C-B at x = —d/2 and B-TZ at x = d/2 (figure[2l). 

While the Schrodinger equation determines the functional form of the wavefunctions 
in the three regions, their actual shape has to be determined by the boundary conditions. 
Since tunneling is viewed as scattering of an incoming electron at the barrier (figure P), 
scattering boundary conditions have to be applied. Note that only electrons in occupied 
states of the source electrode can tunnel into unoccupied states of the drain electrode, due 
to the Pauli exclusion principle. 

Assume for the following that the wavefunctions in the leads and in the barrier at 
energy E have been computed according to the approaches introduced in section HI Since 
E is fixed, the explicit energy dependence will be dropped. 



5.1 Nonmagnetic leads 

5.1.1 Determination of the total wavefunction. In the present one-dimensional 
model, there is only a single Bloch state in each lead with positive (or negative) velocity. 
Thus, the total wavefunction in £ is superposed by the incoming Bloch state ipi{x) (with 
velocity Vi > 0) and the refiected Bloch state rx(jj.{x) (with velocity Vj- = — fj), V'i(2^)+^'^r(a;). 
In 71 there is only the transmitted Bloch state tipt{x) with f t > (figure [l]). Within the 
barrier B, two exponential functions are taken, Aexp^Kx) + B exp{—Kx). The matching 
conditions provide the probabilities for transmission through and refiection at the barrier 
given by |tp and |rp, respectively. 

Within Landauer-Biittiker theory, reservoirs contribute each incoming scattering chan- 
nel with the same current ji. Therefore, either the Bloch states [12] or the transmission 
probabilities have to be normalized to unit current, 

T=\tf-, (11) 

Ji 

with the current jt carried by being calculated from j = lm{-ip* -^"tp)- The transmission 
T measures the transmitted current which is carried by ipt upon feeding with ipi. Current 
conservation implies (1 — |rp)ji = |tpjt- 



5.1.2 Current and conductance. Applying a bias voltage l^ias to the tunnel junction 
shifts the Fermi energies of the electrodes, as is shown schematically in figure El For zero 
bias, the Fermi energies Ep of the electrodes C and TZ coincide (E^ = E^, figure Eti). For 
Vbias > 0, the potential V^{x) is shifted rigidly to lower energies, which is achieved by 
taking as potential average Vq = —Vbias instead of Vq = (as being fixed for £). Evidently, 
both the band structure and the Fermi energy of 7^ are shifted alike, E^{k) = E^{k) — Vua.s 
and Ep- = Ep — Vbias (figure Ob). Hence, electrons can tunnel from £ to 71. Since 
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only electrons of occupied states in L (with energy E < Ep ) can tunnel into unoccupied 
electronic states in TZ{Ep — Vbias < E), an energy window opens up in which tunneling can 
take place {E e [Ep — Vbias,-£'F]) cf. the orange areas in figure IH]). 

The total current / through the tunnel junction, which flows from the source electrode 
£ to the drain electrode TZ, can be expressed for small Vbias (linear response) as / = G Vuas. 
Within Landauer-Biittiker theory, it is given by the sum over all transmission probabilities 
T{E) in the energy window [Ep — Vbias, Ep], 



In general, T{E) implies a sum over the transmission probabilities of all incident Bloch 
states ipi and all transmitted Bloch states ipt- In the present one-dimensional model, with a 
single electronic state available in each electrode, T{E) reduces to a single term [cf equation 
(fTI|) ]. The quantum of conductance Gq is in in atomic units (2e^//i), where the factor 2 
comes from the spin degeneracy. In experiments, the differential conductance d//dVbias is 
typically recorded, which is computed numerically from equation (fT2jl . 
The conductance G of the tunnel junction. 



depends on details of the tunnel junction via the transmission T{E), e. g., on the electrode 
band structures and properties of the tunnel barrier. 

5.2 Ferromagnetic leads 

In nonmagnetic systems, each electronic state can be occupied by two electrons, one with 
spin up, the other with spin down. The spin-dependent band structures of each lead 
coincide [E^{k) = E{k) and E^{k) = E{k) — Vbias] and the number of spin-up electrons 
equals the number of spin-down electrons. This degeneracy is lifted in magnetic systems, 
and the spin-dependent band structures are rigidly shifted by the exchange splitting Vex > 
with respect to each other, i. e., E^{k) = E{k) =f Kx/2 and E'^{k) = E{k) - \4ias =F Vex/2 
(13)). Since the electronic states are occupied up to the Fermi energy E-p, there are usually 
more spin-up electrons than spin-down electrons, which leads to the notation majority 
(minority) electrons indicated by "+" ("—")■ 

Associating majority (minority) with spin-up (spin-down) electrons in C [21], defines 
a net magnetic moment M"^ which is parallel to the magnetic moments of the electrons, 
due to the surplus of majority electrons. In a magnetic tunnel junction two configurations 
remain for the right lead. Either the net magnetic moment in TZ is parallel to that 
in C (M"^ II M^) which is called P configuration or anti-parallel (M^ || — M''^) which is 
called AP configuration, respectively. 

In experiment, the configurations can be switched by an external mag Because the spin 
is conserved during the tunneling process, the total current in each configuration is given 




(12) 




(13) 
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by the sum of the two spin-resolved currents ("two current model"). In P configuration, 
with electrons tunneling from majority (minority) states into majority (minority) states, 
one has the upper situation of figure [SI 

Jp = /++ + /__, (14) 

and accordingly for the AP configuration, with majority and minority interchanged in TZ 
(bottom in figure EJ, 

Iap = I+- + /-+. (15) 

The four spin-dependent currents are related to the conductance contributions Gp = 6*++ + 

G and Gap = G^^ + The quantum of conductance Gq is now 2tt (i. e., /h). For 

a nonmagnetic junction, all spin-resolved currents are identical (/++ = / = /_+ = /-+). 

In typical tunnel junctions in P configuration, one of the spin-resolved currents exceeds 

the other by far (for instance, /++ ^ I in figure [HI), whereas in AP configuration, both 

currents are of almost the same size (/+_ ~ The tunnel magneto-resistance (TMR) 

5 is a measure for the spin-dependent transport and defined here by the asymmetry of /p 
and /ap, 

c- -^p — -^AP Cp — Gap „N 
^= T = r ■ ^^^^ 

Jp + J-AP <-^P + L^AP 

The tunneling calculations for ferromagnetic leads proceed as those for nonmagnetic 
leads but with the band structures shifted by =fKx/2, giving rise to the four currents I±±. 
Eventually, Jp, Jap, and 6 are computed from the latter. 

6. Results and discussion 

The theory of tunneling, as described in the preceding sections, was implemented in a set 
of Mathematica™ notebooks [8]. Because the results for a magnetic tunnel junction 
(subsection I6.2|l can easily be understood from those for a nonmagnetic one, the latter are 
presented and discussed first. For the following, the basis set is Q = {—^g, —'2.g, . . . , +2g} 
(r-Q = 6). The Fourier coefficients Vg of the lead potentials are chosen as Vq = OH, 
Vg = 0.012 H, V2g = 0.008 H, and Vsg = 0.004 H. 

6.1 Nonmagnetic leads 

6.1.1 Nearly- free vs free electrons. The Hamilton matrix for free electrons is diag- 
onal {Vg = 0), and the Bloch states become pure plane waves, (f)^'^^ = exp[i(fc + q)x], with 
band structure E^'^\k) = {k + g)^/2 and velocities v^'^\k) = k + q. The density of states 
decreases with energy, D{E) = 27idk{E)/dE = 7c^/2/E. 

The nonzero potential in the NFE case results in band gaps at = and k = ±g/2. The 
band gaps at k = ±g/2 are between band 1 (lowest in energy, black in figure andEt) and 
band 2 (orange) and show up between 0.112 H and 0.137 H. The first gap at A; = appears at 
0.492 H between bands 2 and 3. The band widths are much larger than the band gaps. For 
band 1, for instance, the band width is EW(^/2) - ^^(0) = 0.112 H + 0.006 H = 0.118 H, 
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whereas the band gap is E^'^\g/2) — E^^\g/2) = 0.025 H. Thus, the electrons behave as 
nearly free (NF), that is, the electronic properties differ from those of free electrons in the 
small energy ranges at the band gaps. 

The velocities v^^\k) = dE^^\k)/dk for NF electrons (figure IHlb) which are linear in 
a large part of the BZ, vanish at the band edges {k = ±g/2 and k = 0). Further, they 
change sign at k = 0. Since electrons with zero velocity cannot contribute to the current, 
a considerable effect on the tunnel current is expected at the band gaps. The DOS D{E) 
(figure Et) also resembles that of free electrons but becomes singular at the band edges. 
These Van-Hove singularities will also show up in the tunnel current. Considering the NFE 
Bloch states, the expansion coefficient Cq^'* of band 1 dominates in a large part of the BZ, 
that is, the wavefunction is to a very good approximation equal to the single plane wave 
with q = (figure int). At k = ±g/2, other plane waves become mixed in: \clI\ > at 

k = +g/2 and |cg^''| > at k = —g/2. Band 2 shows the same qualitative behaviour, but in 
addition mixing at k = 0. In summary, differences between free and nearly- free electrons 
in the leads show up close to band gaps. 

In order to work out how the tunneling is modified with respect to that of free electrons, 
the transmissions T{E) vs barrier width d are compared for energies close to and well 
below the band gap at k = g/2. For d = 0, the transmission is perfect in all cases {T{E) = 
1] because the incoming Bloch states are not refiected at the barrier (figure Et^)- With 
increasing d, T{E) decreases due to the decaying wavefunctions in B. For E = 0.112 H, 
i.e., close to the band edge, the transmission for free electrons is significantly larger than 
for NF electrons, as can be explained by their larger velocity. On the contrary, there is 
no apparent difference for E well below the band edge (here: E = 0.050 H) because in 
both cases the Bloch states comprise a single plane wave. In summary, band gaps have 
a pronounced effect on the transmission and show also up in the bias dependence of the 
current. 

6.1.2 Bias dependence of tunneling. Before turning to current and conductance, 
the dependence of the transmission on bias voltage and barrier width is addressed for 
fixed energy (figure 03) . The overall smooth shape is disturbed by sharp minima at about 
0.013 H and 0.4 H bias which are related to the band gaps in TZ. The monotonous decay 
with barrier width d is consistent with that in figure [3 

At first glance, one would expect that a current / through a tunnel junction increases 
with bias voltage, in accordance with Ohm's law. This behaviour can indeed be observed 
in figure IHti, but only as a general trend. For bias below 0.1 H and at about 0.4 H ample 
deviations occur, with I showing even a decrease. These features show up more pronounced 
in the differential conductance dl/dV (figure |Hlb), where the linear dependence of I on bias 
can be clearly resolved between 0.10 H and 0.35 H. The conductance drops rapidly with 
bias but remains almost constant for Vbias > 0.1 H (figure Et). 

The aforementioned features can be discussed in terms of band structures and transmis- 
sion probabilities. However, a simpler access is provided by the number of states within the 
energy window of tunneling [Ep — Vbias, E^], an approach to be viewed as an extension of 
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the Maekawa-Gafvert model to nonzero bias (section [2)) . The number A^^(Vbias) of incident 
states in L is 



where the factor 1/2 takes into account that only half of the Bloch states show positive 
velocity. For within band 1 (at 0.1 H), N-^ increases monotonously until Vbias = -Ep — 
-E^ot, i-e., when the bottom of band 1 is reached [-Ebot = -E^^-'(O); figure |Hll]. Because 
A^^ remains constant for larger bias (figure El and figure Ej), it is sufficient to restrict the 
integration in equation (fTTjl to [-Emin^ -^f]) with £'^in = i^cLx(i?bot5 -^f ~ Hias)- 

The number iV/^(Vbias) of transmitted states in 7?. depends on the energy range in which 
incident electrons are provided by L. Thus, 



As is evident from figure [Hll, the general shape of Nl^ consists of an increase for bias up 
to about 0.13 H and a slight decrease for larger bias. But on top of this, it shows a richer 
structure than . Upon increasing Vbias, a kink appears either when a band edge enters 
the energy windows (at i?p) or when it leaves the window (at £'min)- This is supported by 
considering the DOS of Tl. at the boundaries of the energy window from (figure |Hfe). At 
the upper boundary [D^(i?p), black in figure [Ht], a band gap produces a decrease of the 
current since there are less transmitted states available. The lower boundary \D^{E^^^^ 
orange], however, shows much less effect. The small kink in / and G at 0.1 H is due to A^j^, 
because it appears just when — E^^^ = Vbias- 

In summary, the structures in the bias dependence of both the current and the conduc- 
tance can be explained by the density of states at the boundaries of the tunnel-energy win- 
dow and by the number of incoming and transmitted states. The actual current, however, 
is further influenced by the tunneling probability which appears to be of minor importance 
due to its smooth shape. 

6.2 Ferromagnetic leads 

Having investigated the transport properties of tunnel junctions with nonmagnetic leads, 
we now turn to junctions with ferromagnetic leads. With regard to subsection 15. 2| the 
spin-dependent currents I±± and conductances G±± were obtained from calculations for 
nonmagnetic leads in which the band structures E±{k) were rigidly shifted by the exchange 
splitting =fKx/2 (chosen as 0.03 H in the following). 

Considering the conductances for the P configuration (figure [Ht), the arbitrarily chosen 

exchange splitting causes a bias dependence of the conductance contribution G which 

is very similar to the nonmagnetic case (figure [Hb) . The conductance contribution 6*++, 
however, is zero for small bias until the bias bridges the gap to the top of band 1 in £. For 

larger bias, G++ is finite but always much smaller than G Consequently, Gp is mainly 

determined by G (figure [Hti) ■ 




(17) 




(18) 
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In comparison to the P configuration, the conductances in AP configuration appear 
slightly modified. G-+ is zero until Vuas bridges the gap to the bottom of band 2 in TZ. 
It increases until a maximum is reached, drops rapidly, and remains almost constant for 
Vbias > 0.1 H until the next gap occurs. (^-^-(AP) is similar to G (P) but the correspond- 
ing extrema are shifted to lower energies by Vex- G+-(AP) is similar to G++(P) with the 
same onset at bias voltages reaching the top of band 1 in £. Gap is determined by 

The tunnel magneto-resistance 5 (equation fT6|) shows a pronounced bias dependence 
which is related to the bandstructure, especially to the position of the bandgaps of the leads 
(figure Et). First of all, a zero-bias anomaly shows up which is caused by the half-metallic 
behaviour of the leads. This is — in contrast to other explanations [22] — a pure electronic 
effect. In general, 5 decays with increasing Vbias except the infiuence at the bandgaps. As 
discussed above (figure ^ , a bandgap in TZ at the upper boundary of the energy window 
causes a decrease of the conductance. Minima in Gp manifest themselves as minima in 6. 
Conductance minima at Gap (G_+) cause maxima at Vbias = 0.036 H and 0.39 H. These 
maxima and the subsequent minima are separated by Kx- 

Concerning the thickness dependence of the tunneling magneto-resistance: the absolute 
values of 5 increase. Hence, the characteristic features of 5 become more pronounced. The 
reason for this behavior is related to the fact that Gp + Gap drops faster with increasing 
barrier thickness than Gp — Gap- 

7. Concluding remarks 

The main purpose of this paper was to demonstrate the role of a realistic bandstructure of 
the leads, described within a NFE model, for the current-voltage characteristics of tunnel 
junctions. The consideration was extended to ferromagnetic junctions in order to discuss 
the bias dependence of the tunnel magneto-resistance. In difference to a free-electron 
model a pronounced TMR effect was obtained. The TMR reveals a zero-bias anomaly of 
electronic origin caused by the half-metallic behaviour of the leads. Driven by the position 
of the band gaps in the drain electrode, with respect to the Fermi energy in the source 
electrode, pronounced positive and negative TMR values occur. 

As further extensions to the presented theory, one could conceive to discuss the spin 
polarizations defined in terms of the number of electrons in the leads and of the spin- 
resolved currents I±±. Contact to the popular Julliere and the Maekawa-Gafvert models 
[9, 10] can be made by trying to explain the TMR 5 by these spin polarizations. However, 
one has to be aware that a possible agreement does not hold in general for three-dimensional 
models because in one dimension the conductance is related simply to the number of states 
in the leads. Another issue could involve different leads which result in a nonsymmetric 
current- voltage characteristics [/(Vbias) 7^ —/(—Vbias)]- 
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Figure captions 



Figure 1: Sketch of a tunnel junction with ferromagnetic leads (£, TZ) and a step-shaped 
barrier (B). An electron incident from the left (with current ji) is reflected (current jr) 
and transmitted through B into TZ (jt). In the leads, the potential depends on the spin of 
the electron: for a spin-up electron it is lower (t, dashed red line) than for a spin-down 
electron (J,, solid red line). 



Figure 2: Potential V{x) of the tunnel junction. For the barrier B, a step-shaped potential 
with height Vban- and width d is assumed. The potential in the leads C and 71 is periodic 
with the lattice spacing a. The interfaces between the three regions are at a; = ±d/2. 



Figure 3: Sketches of the electronic structure of a tunnel junction with a bias voltage 
applied, as depicted schematically by the density of states (black areas) for the leads JC 
and 71. For zero bias, the Fermi energies Ep coincide (a), whereas for positive bias Ep- is 
lowered to E^ — Vu^s in T^, opening up the energy window of tunneling [orange area in 
(b-d)]. The size of the energy window increases with T4ias (b) until its lower boundary 
Ep — Vuas reaches the bottom E^^^ of the >C-bands (c, d). The barrier is denoted B. 



Figure 4: Spin-resolved band structures in the ferromagnetic leads of a tunnel junction in P 
configuration. The majority-electron bands (+, black) are shifted by > with respect 
to the minority bands (— , orange) in both leads. The junction is biased by Vuas > 0. 

Figure 5: Configurations of a tunnel junction with ferromagnetic leads. The lead magneti- 
zations are either parallel (M"^ || M''^; P configuration, top) or anti-parallel (M"^ || — M^; 
AP configuration, bottom). The spin-resolved currents I±± are visualized by horizontal 
arrows. 



Figure 6: Electronic structure of a nonmagnetic lead, with a focus either on the E(k) (left: 
a-d) or on the k{E) description (right: e-h). (a) Bandstructure E{k) of the lowest (black, 
band 1) and the second (orange, band 2) band in the first Brillouin zone, k e] — g/2, g/2]. 
(b) Velocity v{k) of the corresponding Bloch states, (c, d) Fourier coefficients Cq\k) of 
the Bloch states associated with band 1 (c) and 2 (d). (e) Bandstructure, but in k{E) 
representation, (f) Density of states D{E) for the bands shown in (a) and (d). (g, h) 
Fourier coefficients of (c) and (d), but for /c > (g) and A; < (h). 



Figure 7: (a) Transmission T through an unbiased tunnel junction with nonmagnetic leads 
(Harr = 0.12 H). For energies close to {E = 0.112 H; black lines) or well below {E = 0.050 H; 
orange lines) the band gap at k = g/2 (E = 0.112 H), T{d) is plotted vs barrier width d 
for nearly-free (solid) and free (dotted) lead-electrons, (b) Transmission probability T{E) 
vs bias and barrier width d at E — 0.1 H. 



Figure 8: Current-voltage characteristic of a tunnel junction with barrier width d = 2.5 a 
and nonmagnetic leads, (a-c) Total current /, differential conductance d//dVbias, and 
conductance G vs bias voltage, (d) Number of states incoming from £ (black) and outgoing 
in TZ (orange) . (e) Density of states of the right lead 71 at E^ (black) and at E^^^ (in units 
of the lattice constant). The Fermi energy is 0.1 H. 



Figure 9: Spin-dependent tunneling for Ep — 0.1 H and exchange splitting Vex — 0.03 H. 

(a) Spin-dependent conductances G++ and G for the P configuration vs bias for barrier 

width d = 2.5 a. Their sum is denoted Gp. (b) As for (a) but for the AP configuration, 
with Gap = G+^ + (c) Tunnel magneto-resistance S obtained from the conductances 

of (a) and (b). (d) Dependence of the TMR 5 on barrier width, as indicated by different 
line styles. 




Figure [T] 




Figure [2l 



c B n 



^hot 



(a) zero bias 



W E'^-Fbias>Ett 



c B n 



— Vbias 



a 



bot 



c B n 



Ep — Vbias 



(d) E^ - Vbias < E 



bot 



Figure El 




Figure [HI 




Figure El 




Figure 



Figure 133 








0.1 



0.2 



0.3 



0.4 



0.5 



o 



0.02 









: 1 


' 1 ' 1 ' : 


1 _ 








1^ (b)- 


: 1 


, , , , , ^ 


A. 




^bias ^■^-^ 



(c) 



(d) 



(e)- 



0.5 



Figure |H1 



0.003 




oo 



1 

0.5 



-0.5 



1 1 1 


1 ' 1 ' 1 ' 

— d=la 

— d = 2 a :\ — 
■ ■ ■ ■ d = 3 a 
















'•>■ 


— m w" 




1 1 1 


1,1,1',- 







0.1 



0.2 



0.3 



0.4 



0.5 



Figure IHl 



